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ABSTRACT 

The acceleration of an outflow along inclined magnetic field lines emanating from an 
accretion disc can be studied in the local approximation, as employed in the compu- 
tational model known as the shearing box. By including the slow magnetosonic point 
within the computational domain, the rate of mass loss in the outflow can be calcu- 
lated. The accretion rates of mass and magnetic flux can also be determined, although 
some effects of cylindrical geometry are omitted. We formulate a simple model for the 
study of this problem and present the results of one-dimensional numerical simulations 
and supporting calculations. Quasi-steady solutions are obtained for relatively strong 
poloidal magnetic flelds for which the magnetorotational instability is suppressed. In 
this regime the rate of mass loss decreases extremely rapidly with increasing field 
strength, or with decreasing surface density or temperature. If the poloidal magnetic 
field in an accretion disc can locally achieve an appropriate strength and inclination 
then a rapid burst of ejection may occur. For weaker fields it may be possible to study 
the launching process in parallel with the magnetorotational instability, but this will 
require three-dimensional simulations. 
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1 INTRODUCTION 



In a well known paper. iBlandford &: Pavnd ([1983) presented 
a mechanism by which jets or winds can be launched from 
accretion discs. Consider a disc that is threaded by a poloidal 
magnetic field and is sufficiently ionized for ideal magneto- 
hydrodynamics (MHD) to be a good approximation. Above 
the surface of the disc, where the magnetic field is dynam- 
ically dominant, gas tends to rotate with the same angular 
velocity as the part of the disc to which is magnetically con- 
nected; its motion in the meridional plane also tends to be 
parallel to the magnetic field. Above a thin Keplerian disc, 
the net acceleration of gas parallel to the field is directed 
away from the disc if the poloidal field is inclined at more 
than 30° to the vertical, when both the centrifugal and grav- 
itational forces are taken into account. 



accretion, accretion discs - ISM: jets and outflows - magnetic flelds 



that are directed away from the rotation axis are suitable for 
launching jets. 



2 LOCAL APPROXIMATION 
2.1 Introduction 



Fig. [T] based on Fig. 1 of IBlandford fc Pavnd (|l982h . 
illustrates the effective potential experienced by gas that is 
forced to corotate with the angular velocity of a Keplerian 
orbit of radius ro in the plane z = Q. Near the saddle point 
of the effective potential at this same location, acceleration 
away from the disc can occur if the inclination exceeds 30° 
either to the left or to the right. Indeed, the picture is sym- 
metrical near the saddle point. Globally, however, only flows 



The shearing sheet (Fig. [2 iGoldreich fc Lvnden-B"eiilll965l 'l 
is a widely adopted local model of an astrophysical disc. The 
sheet is centred on a reference point that follows a circular 
orbit around the central mass, and the frame rotates with 
the angular velocity Q. of this orbit. The differential rota- 
tion of the disc is represented locally as a uniform parallel 
shear flow in the rotating frame. The local approximation 
has two symmetries not present in the global system: it is 
invariant under translation in the x direction (allowing for 
a Galilean transformation in the y direction to correct for 
the velocity shift) and under rotation by 180° about the z 
axis. The latter symmetry might appear to make the shear- 
ing sheet unsuitable for studying jet launching, because it 
fails to distinguish between radially inward and outward di- 
rections; however, as note d above, this symme try is indeed 
present in the problem of IBlandford fc Pavng (|l98^ ) close 
to the disc. 
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Figure 1. Contours of the effective potential experienced by mat- 
tor tliat is forced to rotate with the Keplerian angular velocity at 
radius tq. The contour values are unequally spaced. Dotted con- 
tours correspond to values lower than that of the saddle point. 
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Figure 2. The shearing sheet, or local approximation. Coordi- 
nates X, y and z are measured in the radial, azimuthal and vertical 
directions. 



For simplicity, we assume in this paper that the disc 
orbits in a Newtonian potential, giving rise to Keplerian 
rotation in the absence of other radial forces, and that the 
gas is isothermal, its pressure and density being related by 
p = c^p, with a uniform isothermal sound speed Cs. The gas 
satisfies the equation of motion, 

pl-^+u- Vu + 20 e^ X IX j = -pV$ - Vp 
+J X B + V T, 

the equation of mass conservation. 



(1) 



| + V.(pn) = 0, 



the induction equation 

dB 

~dt 



V X {ux B - riV X B) 



(2) 



(3) 



and the solenoidal condition 

V ■ B = 0. (4) 

Here $ = ^fl^{z^ — 3x^) is the effective (tidal) potential 
in the local approximation, which comes from expanding 
the sum of the gravitational and centrifugal potentials to 
second order in the distance from the centre of the sheet. 
Also J = fig^'V X B is the electric current density and Tij — 
piy{uij +Uj,i) + p{ui, — ^i')uk,kSij is the viscous stress tensor. 
We assume that the kinematic shear and bulk viscosities v 
and fb and the magnetic diffusivity r/ are uniform. 

The hyperbolic contours of $ in the xz plane agree with 
those of the effective potential plotted in Fig. [T] close to 
the saddle point and reproduce the critical inclination of 
30°. The reason for this is of course that $ is the effective 
potential experienced by matter that is forced to rotate with 
the Keplerian angular velocity Q, of the reference point. 

The basic state of the shearing sheet in the absence 
of magnetic fields consists of the Keplerian shear fiow u — 
~^Qxey together with the hydrostatic density distribution 
p = Po exp{—z^/2H^), where po is a constant and H = Cg/fi 
is the isothermal scaleheight. The departure from the basic 
Keplerian fiow is denoted hy v — u + ^^x Sy. 

This system of equations can be solved numerically in a 
finite shearing box (e.g. iHawley, Gammie fc Balbus 1995 ), 
in which case shearing-periodic horizontal boundary con- 
ditions apply to the solution (to v rather than u), while 
various vertical boundary conditions are permissible. The 
shearing box has been widely employed in treatments of 
the m agnetorotational instability (MRI) l|Balbus fc HawlevI 
Il998l . and references therein), but has not been used for stud- 
ies of jet la unching. [The recent work of I Suzuki fc Inutsukal 
(2009) and lSuzuki. Muto fc Inutsukal (|2010l ') finds mass loss 
from the computational domain but does not consider the 
systematically inclined fie lds relevant for the mechanism of 
iBlandford fc Pavnd l|l982h .] 



2.2 Vertical boundary conditions 

For any choice of vertical boundary conditions at z — ±Z, 
and with shearing-periodic horizontal boundary conditions, 
the horizontal average over the box of B^ is independent of 
z and t, and is determined by the initial conditions. Within 
the local approximation, the type of poloidal magnetic field 
configuration that is favourable for jet launching consists of 
a uniform vertical field Bz together with a radial field Bx 
that is odd in z and tends to a non-zero constant at large z; 
this represents a field that is straight, inclined, current-free 
and therefore force- free in the low-density gas at large \z\ 
and bends symmetrically as it passes through the disc, pro- 
ducing a radial Lorentz force through the azimuthal current 

Jy (Fig.©. 

In this type of study, where a thin disc is to be connected 
to a jet, the question arises of which quan tities are deter- 
mined by the disc and which by the jet fsee lOgilvie fc Livid 
|2001, and references therein). Efficient outfiows pass through 
a slow magnetosonic point not far above the surface of the 
disc and through an Alfven point much higher up (|Spruitl 
Il996l ). It is convenient to think of matching the disc to the 
jet a small distance above the slow magnetosonic point, in 
a region where the magnetic field is predominantly poloidal 
and approximately force-free. Local studies of the disc can 
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dt I (pvy) dz + -O, / (pvx) dz 
J-z ^ J -z 



Figure 3. Geometry of the poloidal magnetic field in the (a;, 2) 
plane that is favourable for jet launching in the local approxima- 
tion. The field lines could equally well bend the other way. 



focus on the dynamics of the disc and the passage of the 
outflow through the slow point, which determines the rate of 
mass loss. However, the rate of angular momentum loss, and 
therefore the magnetic torque on the disc, are determined by 
the passage of the outflow through the Alfven point, which 
usually lies outside the domain of the local approximation. 
Since the magnetic torque is proportional to the azimuthal 
magnetic field, the value of By in the matching region is de- 
termined by the jet region, not the disc region. Given the 
distribution of poloidal magnetic flux on the midplane of 
the disc, the global shape of the poloidal magnetic fleld is 
also determined by the jet region. It is therefore appropriate 
to impose on the local disc the value of B^, and the values 
of Bx and By in the matching region. This can be done by 
a novel type of boundary condition in which the horizontal 
components of the magnetic field are specified at the vertical 
boundaries, with equal and opposite values at the top and 
bottom because of the desired symmetry. At the same time, 
the vertical boundary conditions should leave the density 
and the velocity free to evolve according to the dynamics of 
the outflow that develops. 



ByBz 

pVyVz + Ty 

MO 



(6) 



In most shearing-box simulations the right-hand sides of 
these equations are negligible or zero, with the result that 
the mean horizontal momentum of the box executes an 
epicyclic oscillation of constant amplitude, which can be 
made to vanish by a suitable choice of initial conditions. 
Here, however, we wish to impose non- vanishing time- 
independent magnetic stresses at the vertical boundaries. 
These provide source terms for the epicyclic oscillation. If 
the initial conditions are chosen carefully, a non-oscillatory 
solution is possible, in which the xz stress is matched by a 
mean flow in the y direction and the yz stress by a mean 
flow in the x direction. The first represents a departure from 
Keplerian rotation induced by the radial Lorentz force of a 
poloidal magnetic field that bends through the disc, while 
the second represents an accretion flow driven by an imposed 
magnetic torque. For most initial conditions, however, a free 
epicyclic oscillation of constant amplitude and phase will be 
superimposed on this non-oscillatory solution. (In fact, de- 
pending on the way that mass loss from a finite computa- 
tional domain is treated, the amplitude may not be constant, 
as is discussed below.) 



3 HORIZONTALLY UNIFORM SOLUTIONS 

We do not attempt the numerical solution of this system 
in multiple dimensions in this paper. Indeed, simulations of 
the MRI including both vertical gravity and a net vertical 
magnetic fi eld present serious nu merical and interpretative 
difficulties (jMiller fc Stonell200Cl : S. Fromang, private com- 
munication; G. Lesur, private communication; J. M. Stone, 
private communication). For the remainder of this paper we 
focus on horizontally uniform solutions that depend only 
on z and t. These are admissible because of the horizontal 
translational symmetry of the local approximation. Then we 
have 






' dVy dVy 1 ^ \ T r, 3 / dVy 



(7) 
(8) 



2.3 Evolution of the horizontal momentum 

Let angle brackets denote a horizontal average over the box. 
Then, by rewriting the horizontal components of the equa- 
tion of motion in conservative form and integrating over the 
box, we obtain 



/Z pZ 

(pvx) dz -2Q (pVy) 
-z J-z 

BxB: 



dz 



A*o 



- PVxVz + T^ 



(5) 



dvz 



dvz 



^l^+"^9. 



\Jx-tjy JyJDx 



-pit z- — 



d 



I , 4 x^W^ 



dp d , . ^ 



as. dEy 



dt dz 



dBy 3 dE^ 



(9) 

(10) 

(11) 
(12) 
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dt 

dB, 
dz 



with 



Ex = v^By — VyBz + rjJx 



Ey — VxBz 



fJ-oJx 



VzBx ~\-r}Jy, 

dBy 
dz 



fJ-oJy = 



dBx 
dz 



(13) 
(14) 

(15) 
(16) 

(17) 
(18) 



Clearly B^ is a constant, which becomes a parameter of the 
problem. The solution of these equations also contains a 
free epicyclic oscillation of arbitrary amplitude and phase. 
That is, the equations are invariant under the transforma- 
tion Vx I— >■ Vx +2asm[Q{t~T)], Vy i— >■ Vy + acos[f7(t — r)], for 
any real numbers a and rljin general, the initial conditions 
will excite this motion at some amplitude. 

In the limit of ideal MHD {v = i/^, = rj = 0) the equa- 
tions form a hyperbolic system, with wave speeds v given 

by 



[(«- 



-A[{i'-i'z)* -{c^+vl){v-Vzf +c^vl 



0, (19) 



where Va ~ (fiop)^^ B is the Alfven velocity. The six solu- 
tions are, of course, the Alfven waves and the fast and slow 
magnetoacoustic waves. 

If the upper boundary z = Z is located above the slow 
magnetosonic point but below the Alfven point, then two of 
the six waves are directed into the computational domain 
and the remaining four waves are directed outward. There- 
fore two boundary conditions should be applied to the six 
dependent variables («; 



X, Vy, Uz, 



p,Bx,By) there, and we do 
this by specifying the values of Bx and By, as discussed 
above. 



4 WAVE SPEED LIMITER 

A difficulty arises in solving the above system of equations 
if the density becomes very small above the disc. The very 
large Alfven speed means that a very short timestep is re- 
quired for stability unless an implicit integration scheme is 
adopted. It is possible to limit the wave speeds in a way that, 
to some extent, mimics the effects of Einsteinian relativity. 
This can be done by introducing a displacement current, 
replacing the definitions of Jx and Jy by 



1 dEx 



dBy 



PoJx + 



t^oJy^^-^ - -1^ 



C2 dt 


d. 


1 dEy 


dB, 


C2 dt 


dz 



(20) 



(21) 



The electric field E undergoes a non-trivial transformation. In 
fact E is not the total electric field because it does not include the 
contribution (proportional to x) from induction by the Keplerian 
shear fiow. 



where c mimics the speed of light, but is a free parameter of 
the nrodel. (Optionally, Jz and Ez can be included as well.) 
When V = Vh = T] = Q we now have a hyperbolic system 
of equations. It is neither Galilean-invariant nor Lorentz- 
invariant, but this is convenient for our purposes because the 
added terms limit the speed of waves relative to the coordi- 
nate systenr (and the numerical grid) while having no effect 
on steady flows and their critical points. We use the wave 
speed limiter nrainly as a convenient way of reaching steady 
solutions using fewer timesteps; these are then also steady 
solutions of the original problem . [The more complicated 
method of lMiller fc Stond l|200(J l does not appear to have 
this property, although it may have other desirable features.] 
In practice the wave speeds are limited to max(c, Cg -I- |mz|). 
Provided that r) is unifornr, the system is still invari- 
ant under the addition of an epicyclic motion, although the 
transformation law of the electric field is modified. 



5 QUASI-STEADY SOLUTIONS 

5.1 ODE system 

The system does not in fact admit truly steady solutions 
with an outfiow unless a source of mass is included. In this 
case we obtain the system of ordinary differential equations 
(ODEs) 



P ^z 



dvx_ 
dz 

AVy 



on 1 d f BxB^ Avx 

2\lVy I = -;- h py—r- 

az \ /io az 



""■■dT + H^E^— 



pVz' 



Avz 
d7 



+ 



Az 



-pQ, z 

{Bl + Bl 
2/io 



AVy 

po az 



. Avz 



+ p{Ui, + ^u)-^_^ 



^(^"^ 



^ d / R ^ dBx 

= -— VxBz - VzBx -I- rj—— 

az V Az 



(22) 
(23) 

(24) 
(25) 
(26) 



2 Az 



VyBz 






ABz^ 
Az 



= 0, 



^ . (27) 



(28) 



where p[z) is the mass source. 

One possibility is to provide a mass source only on the 
midplane z = 0, so that p oc 5{z). The ODEs can then be 
solved in 2 > 0, with p = 0. This procedure resembles the 
way in which steady accretion discs are treated without ac- 
counting for the (slow) increase in the mass of the central 
object. Alternatively, mass can be replenished at a rate pro- 
portional to the local density, so that p = yi^p, where 7 is a 
small positive constant whose value can be adjusted to ob- 
tain a steady solution. In either case, nrass is added to the 
solution at a rate that balances the outflow. 

The solution is expected to be symmetrical about the 
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midplane, with p, Vx and Vy being even in z while v^, B^ 
and By are odd. In the case of a mass source locahzed on 
the midplane, however, v^ is non-zero (and discontinuous) 
at z — 0. 

The ideal MHD problem (in which u = Vh = r] = Q) has 
critical points wherever the wave speeds v given by equa- 
tion (|19[) vanish. We assume that p > and iiz > for all 
z > 0. We expect that the flow will accelerate beyond the 
slow magnetosonic speed a few scale heights above the mid- 
plane, if the inclination of the poloidal field is appropriate. 
In a successful outflow the Alfven and fast magnetosonic 
speeds will also be exceeded, but this will usually happen at 
a greater distance beyond the range of validity of the local 
approximation. We therefore consider only the slow magne- 
tosonic critical point here. 

We apply the following boundary conditions. At the top 
of the domain, z = Z, representing the 'surface' of the disc, 
we prescribe the values of Bx and By to be B^b and By^. 
The constant Bz is also given as an input parameter. The 
inclination angle i of the poloidal magnetic field is defined 
through Bxa ~ -Bz tan i. At the midplane, z — 0, the sym- 
metry conditions Bx = By = apply, as well as «z = in 
the case of the distributed mass source. Finally, we require 
the surface density 



pdz 



(29) 



to equal a given value, which is one of the parameters of the 
model. 

Of particular interest as an output of the calculation is 
the rate of mass loss per unit area from the upper surface 
of the disc, 



■ pvz 



(30) 



The timescale for depleting the mass of the disc is E/(2mw). 
Radial mass transport can also be induced in this 
model, despite the symmetries of the local approximation. 
Although the xy components of the viscous, Reynolds and 
Maxwell stresses are independent of x and therefore do not 
drive an accretion flow, the magnetic stress ByB^/po act- 
ing at the vertical boundaries can do so. In a steady state, 
the integrated azimuthal momentum equation ([6]) relates the 
radial transport velocity of mass, Vm, defined through 

Sii,„ = / pvx dz, (31) 

J -z 

to the yz stresses acting at z = ±Z. 

Furthermore, the radial transport velocity of poloidal 
magnetic flux, v^, is given by 



Ey = VxBz 



VzBx + r;—- — = v^Bz 
az 



(32) 



The physical significa nce of this quantity is discussed by 
iQgilvie fc LMol (|200ll ) and Gullet & Ogilvie (MNRAS, sub- 
mitted). It contains both an advective contribution, from 
motion in the meridional plane across the field lines, and a 
diff usive contribution. 

lOgilvie fc Livid (|200ll ) solved a closely related problem 
for the vertical structure of magnetized discs with outflows. 
There are several differences. In that work, the problem was 
separated into an optically thick disc and an optically thin 
atmosphere. In the disc, the vertical velocity was neglected 



and the temperature was determined according to an energy 
equation including viscous and resistive heating and radia- 
tive diffusion. In the atmosphere, an isothermal outflow was 
computed along rigid magnetic field lines. Another difference 
is that the vertical transport of momentum by viscosity was 
neglected, while the radial transport of angular momentum 
was modelled in a way that goes beyond the standard local 
approximation. 



5.2 Ideal MHD problem 

We consider here the ideal MHD problem {v — V]^ = rj = 0) 
with mass replenishment at the midplane. Assuming that 
the inclination of the poloidal magnetic field exceeds 30°, 
we expect to find a slow magnetosonic point at some height 

Z — Zs- 

Given that m„ = pv^ and Ey = VxB^ — v^Bx are both 
independent of z, we can choose our vector of dependent 
variables to be 



X 



(33) 



The quantities p and Vx follow from the values of m„ and Ey, 
which have the nature of eigenvalues, while B^ is a constant 
parameter. The remaining ODEs can be cast in the form 

. dX ,, 



dz 



where 



A = 






BxVz 


vi 


-vL 





Vz 










MOP 





Vz-^ 




Sx 


By 








MOP 


MOP 


-if. 


By 







Vz 



(34) 



(35) 



and 



2D,BzVy 

— h^Vx 



'^QBx 



(36) 



At the slow magnetosonic point where 



:(Cs +«a) 



r(^ 



+ ^a) 



1/2 



1/2 



(37) 



the matrix A is singular. For a regular solution that passes 
smoothly through the slow point with flnite derivatives, we 
require LY = 0, where 



L^ Bx 



ByBz 



Bl 



P-OP^z 



ByVz\ 



is the null left eigenvector of A satisfying LA. 
point. The regularity condition is therefore 



(38) 
at this 



ByBzVx-4:BxBzVy+WxByVz+2D.z{B:-p.opv:) = 0. (39) 

A method of solution is as follows. The quantities z^ 
and Ey are guessed, as well as the values of Bx , By and p at 
the slow point. The value of Vz at the slow point follows by 
definition and hence the quantity rh^. The regularity condi- 
tion determines the value of Vy . This is sufficient information 
to integrate the ODEs in each direction away from the slow 
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point. Away from the slow point and other critical points, 
the matrix A can be inverted to find dX /dz from Y . We 
integrate down to z — and require that Bx = By — 



there. We also integrate up to 



Z and require B^ and 



By to equal specified values there. A fifth condition is given 
by the mass integral (I29|l . Using Newton-Raphson iteration, 
the five guessed quantities are adjusted to meet the five con- 
ditions. 

This method fails if the Alfven point is encountered 
below z — Z, because then the matrix A is singular and 
the solution cannot be made to pass smoothly through the 
critical point and match the specified values of B^ and By 
at z — Z. 

A reduced version of the method is to assume that B^ 
and By have their limiting values already at the slow point. 
In that case no upward integration is required. There are 
only three quantities to guess and three conditions to be 
met. 



5.3 Simplified version 

A simplified version of the problem is obtained by neglecting 
Vx, Vz and By and solving the equations 

B. dBx 



2r2i'„ = 

Mop az 



Q^-^'z-^'^ 



p dz 



2fio 



(40) 



(41) 



this limit the condition for the slow point reduces to v, 
Csi'az/t'a, and the regularity condition (|39|) reduces to 

B^Qz 



2Bx 



(45) 



which agrees with what is obtained by assuming an isoro- 
tational configuration and identifying the slow point as the 
location of the density minimum of the static approximation 
as described above. 



5.4 Numerical simulation 

In many ways it is easier to avoid the technical difficulties 
of directly computing steady solutions and instead to evolve 
the system forwards in time until a steady solution is ap- 
proached, if it is stable. We use the wave speed limiter, usu- 
ally with c = lOcs. A branch of steady solutions can be 
followed quasi-statically by slowing changing a parameter 
(such as Bz) during the simulation. 

To integrate the equations we use tenth-order centred fi- 
nite differences in z and a third-order explicit Runge-Kutta 
method in t. The method is usually stable when a version of 
the Courant condition on the timestep is satisfied. Bound- 
ary conditions are implemented with the help of ghost zones. 
In the case of the free 'outfiow' boundary conditions on the 
velocity, the ghost zones are filled using linear extrapolation 
from the computational domain. If desired, mass replenish- 
ment can be carried out by either the localized or the dis- 
tributed method. 



= -^nBx+B.^. 
2 dz 



(42) 



This procedure is equivalent to assuming an isorotational 
configu r ation with a purely poloidal magne tic field, as in 
iQgilviel |l993) and lOgilvie fc LMol (|l998D but with an 
isothermal gas. The boundary conditions are that Bx van- 
ishes at 2 = and has a specified value Bz tan i at large z, 
together with the constraint on the surface density. 
Eliminating Bx, we obtain 



-(-f7. +c,lnp~-vy 







and 



B 



2d^U 



dz 



J- +3Q popVy = 0. 



(43) 



(44) 



< 30° there exist solutions for which p — ^ and 
Bz tani as 2 — > oo, while Vy ~ {^Qtani)z+Vyo, where 
Vyo is a constant. For large z this gives a declining density 



For 
Bx - 



with Inp 



-(1 



)Q^z^/2ct. For i > 30°, the gra- 



dient of Vy is sufficiently large that p would increase with z 
above a certain height, z = [2tan j/(3tan^ j — l)]{—Vyo/fl), 
and eventually Bx would be obliged to vary significantly. 
This indicates a breakdown of the static solution and the 
need for a transcritical outflow. The location of the minimum 
density gives a good indication of the location of the slow 
magnetosonic point in the corresponding solution with an 
outflow, and the rate of mass loss can be est imated from the 
static solution using a correction factor as in lOgilvie fc Livid 
||199^. 

This approximation is expected to be valid when the 
plasma beta at the slow point is small, i.e. v^ ^ Cg. In 



6 NUMERICAL RESULTS 

A dimensionless measure of the vertical magnetic fleld 
strength is 



Bz 



i/poSCsfi 



(46) 



Without loss of generality, we can choose units such that 
/io = E = Cs = f2 = l. The remaining parameters are 
the three magnetic field components {Bxs, Bys, B^) and the 
three diffusivities {v,Vh,ii), whose values in these units can 
be interpreted as Shakura-Sunyaev alpha parameters. 

We focus here on jet-launching solutions with Bxs ~ B^ 
and Bys = 0, for which the inclination of the poloidal mag- 
netic field lines at the upper boundary is i = 45° . The prin- 
cipal remaining parameter is the vertical magnetic field B^. 
(The main effect of imposing a non-zero value of Bys would 
be to apply a magnetic torque to the disc and thereby to 
drive a mean radial fiow.) 

We have computed quasi-steady solutions in ideal MHD 
as described in Section 15.21 using a Runge-Kutta method 
with adaptive stepsize. Mass is replenished at the midplane. 
Solutions are obtained for Bz ^ 0.7; the height of the slow 
magnetosonic point and the rate of mass loss are plotted 
in Fig. |4] In these solutions the inclination of the poloidal 
magnetic field is fixed (i.e. Bx = Bz) at the slow point. 
For Bz > 1 these quantities are accurately predicted by the 
simplified version of the problem described in Sect ion | 5.3I 
This behaviou r is cons istent with that found by lOgilvid 
(1997 ). Ogilvie fc Livid (il998 ) and lOgilvie fc Livid (|200ll ). 
A stronger bending poloidal magnetic fleld produces a larger 
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Figure 4. Top: Height of the slow magnetosonic point, in units of 
Cs/C, as a function of the vertical magnetic field strength, in units 
of ^/^o^CsH, for jet-launching solutions with inclination i = 45°. 
Bottom: One-sided mass-loss rate per unit area, in units of Ef2. 



tializing the simulation with the quasi-steady solution up to 
the slow magnetosonic point. It is straightforward in this 
way to obtain strongly magnetized solutions, e.g. at Bz = 2. 
In this regime the poloidal magnetic field does straighten 
well before 2 = 5, and a steady outflow is achieved, with a 
rate of mass loss that is insensitive to the location of the 
boundary. With a less careful initial condition, a perma- 
nent epicyclic oscillation is generally obtained on top of the 
quasi-steady solution, while other wave modes damp slowly 
through viscosity and resistivity. The epicyclic oscillation 
can be avoided only by initializing the horizontal momentum 
correctly, but its presence has no effect on the jet-launching 
process. 

The presence of non-zero diffusivities slightly reduces 
the rate of mass loss. As an example, for Bz = 2, diffusivities 
of 0.03 reduce the mass loss by about 7% compared to the 
ideal-MHD solution. 

The value of Bz can be decreased slowly and con- 
tinuously during the simulation, or the simulation can be 
restarted periodically with slightly smaller values of Bz ■ As 
this is done, the height of the slow point decreases and the 
rate of mass loss increases dramatically. Eventually the sys- 
tem becomes unstable to the MRI. If the simulation main- 
tains perfect symmetry about the midplane, then the flrst 
MRI mode to be observed is expected to be similarly sym- 
metric; if an antisymmetric seed perturbation is introduced, 
then an antisymmetric MRI mode is expected to set in first. 
In practice it is difficult to observe a clean development of 
the MRI with such an inclined poloidal field. The mass loss 
is so strong in this regime that either the disc empties very 
quickly or the details of the replenishment process become 
important; depending on how momentum is treated in this 
process, a rapid artificial growth of epicyclic motion may 
occur. This transitional regime is also problematic because 
the density of the outflow is so large that the Alfven point 
and fast magnetosonic point approach the disc and the in- 
terpretation of the solutions becomes unclear; there is no 
longer a force-free region of straight fleld lines in which the 
inclination of the fleld can be meaningfully imposed. 

Without mass replenishment, the surface density de- 
creases in time, with the effect that the disc becomes more 
strongly magnetized: the value of Bz/y/fio'^Csfl increases. If 
this parameter starts close to 1, then rapid mass loss oc- 
curs and the parameter increases quickly until the mass loss 
abates. In this way a burst of ejection occurs. In a global 
model of a disc, this condition might be met at different 
times at different radii, leading to a 'wave' of ejection. 



radial Lorentz force, making the disc more sub-Keplerian 
(in the case of outward bending) and presenting a larger 
potential barrier to the outflow. The logarithm of the mass- 
loss rate is proportional to —Bi at large Bz- In practice, 
the outflow is completely suppressed when the dimension- 
less magnetic fleld strength increases by a factor of only 2 
or 3. However, the maximum value of rriw plotted in Fig. [4] 
is about 0.23, which corresponds to the disc being emptied 
on a time-scale of approximately 0.35 orbits! 

We have also carried out numerical simulations in non- 
ideal MHD as described in Section 15.41 In a typical nu- 
merical simulation we solve the equations in the domain 
—5 < 2: < 5 with 1000 grid points and with diffusivities 
u = i/h — ri — 0.03. A clean start can be obtained by ini- 



7 CONCLUSION 

The acceleration of an outflow along inclined magnetic fleld 
lines emanating from an accretion disc can be studied in 
the local approximation, as employed in the shearing-box 
model. By imposing the inclination of the magnetic field at 
the vertical boundaries of the box and resolving the slow 
magnetosonic point within the computational domain, ap- 
propriate solutions with outflows can be obtained and the 
rate of mass loss can be calculated. 

This procedure works best when the magnetic fleld is 
sufficiently strong to suppress the magnetorotational insta- 
bility (MRI). As found in previous work, in this regime the 
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mass-loss rate is extremely sensitive to the dimensionless 
field strength -Bz/VA'oScsfi, allowing the possibility of an 
explosive burst of ejection from any radius in the disc at 
which the appropriate field strength is attained. 

It would be of interest to determine the stability of 
these solutions with respect to perturbations that depend 
on X and/or y, to tes t for the instabilities discussed by 
iLubow fc Spruiti l|l995h and iLubow. Papaloizou fc Pringlj 
1 1994) . For weaker fields it may be possible to study the 
jet launching process in parallel with the MRI. While that 
is a much more demanding problem, the present paper may 
provide some useful guidance. 
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